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The present paper reports high-accuracy cross-section data for the 2 H(n, nnp) reaction in the 
neutron-proton (np) and neutron-neutron (nn) final-state-interaction (FSI) regions at an incident 
mean neutron energy of 13.0 MeV. These data were analyzed with rigorous three-nucleon calculations 
to determine the 1 So np and nn scattering lengths, a np and a nn . Our results are a nn — -18.7 ±0.6 fm 
and o np = -23.5 ±0.8 fm. Since our value for a np obtained from neutron-deuteron (nd) breakup 
agrees with that from free np scattering, we conclude that our investigation of the nn FSI done 
simultaneously and under identical conditions gives the correct value for a nn . Our value for a nn 
is in agreement with that obtained in iv~d measurements but disagrees with values obtained from 
earlier nd breakup studies. 
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The difference in the 1 Sq neutron-neutron (nn) and 
proton-proton (pp) scattering lengths is an explicit mea- 
sure of charge-symmetry breaking (CSB) of the the nu- 
clear force. The high sensitivity of these scattering 
lengths to the nuclear potential strength makes them a 
valuable probe for detecting small potential-energy con- 
tributions like the isospin-dependent forces [[j] [| that 
cause CSB. For realistic potentials, a 1% change in the 
potential strength results in a 30% shift in the scattering 
length. Since it is not technically viable to measure a nn 
directly by conducting free nn scattering experiments, 
one uses few-nucleon reactions that emit two neutrons 
with low relative momentum, i.e., a final-state interac- 
tion (FSI) configuration. The two reactions used to mea- 
sure a nn that have the smallest theoretical uncertainties 
are pion-deuteron capture (7r~ + d — » n + n + 7) and 
neutron-deuteron breakup (n + d — > n+n+p). Curiously, 
measurements using these two reactions give significantly 
different values of a nn : from n~d capture measurements 
the average value for a nn is — 18.6±0.4fm an d from 
kinematically-completc nd breakup experiments the av- 
erage value is -16.7 ±0.5 fin 0. It was suggested that 
the difference in these values has its origin in the three- 
nucleon force (3NF), which would act in the nd breakup 
reaction jjj but not in the other. 

The goal of the present work was to measure a nn 
with the nd breakup reaction to an accuracy better than 
±0.7 fm. To this aim a np was measured simultaneously 
and used as a standard for evaluating our methods and 
for determining the influence of the 3NF in nd breakup. 
We obtained a np and a nn from absolute cross-section 
measurements in the np and nn FSI regions, respectively, 



in nd breakup at an incident mean neutron energy of 
E n = 13.0 MeV. Rigorous nd breakup calculations with 
the Tucson- Melbourne (TM) 3NF potential show that 
the percentage change in np and nn FSI cross sections 
from the addition of a 3NF are equal. Based on this ob- 
servation and the fact that a np has been determined to 
high accuracy by np scattering j^] , we used our extracted 
values of a np to set an upper limit on the 3NF influence 
on the value of a nn determined in the present experi- 
ment. Because 3NF effects might have a stronger energy 
and angle dependence than indicated by our calculations 
with the TM 3NF model, we designed the experiment to 
obtain nn FSI data at the same energy and angle configu- 
rations as that for the np FSI data. Furthermore, guided 
by our calculations of the angle dependence of 3NF effects 
on the NN FSI cross-section enhancement |)| , measure- 
ments were made at several np and nn emission angles. 
The measured cross-section distribution at each np (nn) 
angle was used to determine a value of a np (a nn ), and 
the predicted angle dependence of a np (a nn ) due to 3NF 
effects was investigated. 

All measurements were made using the shielded neu- 
tron source at the Triangle Universities Nuclear Labora- 
tory (TUNL). The experimental setup is shown in Fig. [I]. 
The momentum of the two emitted neutrons and the en- 
ergy of the proton in each breakup event were measured, 
thereby overdetermining the kinematics. The neutrons 
were detected in liquid organic scintillators. The energy 
of the emitted proton was measured in the deuterated 
liquid scintillator CgDi2 (NE-232), which served as the 
deuteron scatterer and will be referred to as the cen- 
tral detector (CD). The energies of the outgoing neutrons 



were determined by measuring their flight times from the 
CD to the neutron detectors. The neutron detectors on 
the right side of the incident beam axis were used for the 
nn FSI measurements. At each angle of the nn pair, one 
neutron was detected in the ring-shaped detector placed 
1.5 m from the CD and the other in the coaxial cylindrical 
detector placed 2.5 m from the CD, which was positioned 
to fill the solid angle of the opening in the ring detector. 




Bicron Detector 

FIG. 1. The experimental setup for the NN FSI 
cross-section measurements in nd breakup. All detectors are 
in the horizontal plane. The two detectors used to measure 
the np FSI cross section at a specific np angle pair lie on lines 
of the same type. The nn FSI cross sections were measured 
at each angle with the pair of ring and cylindrical detectors 
positioned at that angle. 

In the np FSI one neutron moves in the same direction 
as the proton, and the other one is emitted on the op- 
posite side of the incident neutron-beam axis. The neu- 
trons emitted on the right side were detected in either 
the ring-shaped detector or in the cylindrical detector. 
The associated neutrons were emitted to the left side and 
were detected in cylindrical scintillators located as shown 
in Fig. [|. All neutron detectors were filled with a liquid 
scintillator fluid with n — 7 pulse-shape sensitivity (either 
NE213 or BC501A). The active volume of each cylindri- 
cal detector was 12.6 cm dia. x 5.5 cm thick, and that of 
each ring detector was 7.6 cm inner dia. x 13.4 cm outer 
dia. x 4.0 cm thick. The neutron detector efficiencies 
were determined in a dedicated series of measurements 
using neutrons from the 2 H(d, n) 3 He reaction and from a 
252 Cf source. The energy dependence of the relative de- 
tection efficiency of each detector was determined to an 
accuracy of ±1.0%, and the absolute efficiency to ±2.5%. 

The neutron beam was produced using the 2 H(c?, n) 3 He 



reaction. The production target was a 3-cm long cell 
pressurized with 7.8 atm of deuterium gas. The cell was 
bombarded with a 10.0-MeV d.c. deuteron beam, which 
entered the cell through a 6.35-/mi thick Havar contain- 
ment foil and was stopped in a gold end cap. The neu- 
tron energy spread was 400 keV. The detector area was 
shielded from the neutron production target by a 1.7-m 
thick multiple component wall. The neutron beam at 
the CD was defined by a rectangular double-truncated 
collimator, which was designed such that the CD was 
illuminated almost exclusively by unscattered neutrons 
produced in the deuterium gas cell. The deuteron beam 
current on target was about 2 /i^4, and the counting rate 
of the CD was about 400 kHz with a threshold setting 
of one-tenth of the Compton-scattering edge for 7 rays 
from a 137 Cs source. The rate in the CD electronics set 
the limit on the maximum acceptable beam current. The 
pulse-height thresholds on the neutron detectors were set 
at one-third of the 137 Cs Compton-scattering edge. Data 
were collected for a total of 2000 hours. 

The integrated beam-target luminosity was deter- 
mined by measuring the yields for nd elastic scattering 
concurrently with the data from the breakup reaction. 
Since the differential cross section for nd elastic scatter- 
ing can be calculated using realistic NN potentials in the 
Faddeev method with high numerical precision |l0| and 
the calculations agree well with existing data, we elected 
to use calculated cross sections rather than experimen- 
tal data in the luminosity determination. The simultane- 
ous accumulation of the breakup and nd elastic-scattering 
data reduced the sensitivity of our measurements to sys- 
tem deadtimes and absolute detection efficiencies. Also 
this technique eliminated the need for absolute measure- 
ments of the neutron flux and the thickness of the CD. 

The width of the coincidence window in the event- 
trigger circuit was 400 ns and allowed for the concur- 
rent measurements of true breakup events and events due 
to the accidental coincidences between signals from the 
CD and the neutron detectors. All events that satisfied 
conservation of energy within ±2 MeV were projected 
into 0.5-MeV wide bins along the point-geometry kine- 
matic locus. The distance along the kinematic curve of 
E n \ versus E n 2, where E n \ and E n 2 are the energies of 
the two emitted neutrons, will be referred to as S. The 
value of S is set to zero at the point where E n 2 = and 
increases as one moves counterclockwise around the lo- 
cus. The true + accidental and accidental events were pro- 
jected onto the ideal locus separately. The projection was 
done using the minimum-distance technique described by 
Finckh et al. |fl| . The true events were obtained by sub- 
tracting accidental from true + accidental events. 

The measured cross sections were compared to Monte- 
Carlo (MC) simulations that included the energy resolu- 
tion and the finite geometry of the experimental setup. 
The basis of these simulations was theoretical point- 
geometry cross-section libraries generated for a range of 



o-np iflnn) values at incident neutron energies of 12.8, 13.0 
and 13.2 MeV. The point-geometry cross sections were 
obtained from transition matrix elements of the breakup 
operator Uq 

U = (1 + P)f, (1) 

where the T operator sums up all multiple scattering con- 
tributions through the three-nucleon (3N) Faddeev inte- 
gral equation 

f\(j> > = tP\<f> > +(1 + tG )V 4 {1) (l + P)\<f> > +tPG f\<t> > 
+ (l + tG a )v} 1 \l + P)G f\ ( f>> . (2) 

Here Go is the free 3N propagator, t is the AW t-matrix, 
and operator P is the sum of a cyclical and anticycli- 
cal permutation of three nucleons. In the generation of 
the libraries, the terms containing V4, the 3NF potential, 
were set to zero. 

The cross-section libraries were obtained using the 
Bonn-B (OBEPQ) AW potential pf. This potential is 
fitted in the 1 Sq state to the experimentally-determined 
value of a np . The charge-indepence breaking in the So 
AW force is imposed by using for the 1 S'o nn force the 
version of the Bonn-B potential jl2j that was fitted to pp 
data. To account for charge-symmetry breaking in the 
calculations, the total 3N isospin T=3/2 admixture has 
been included |l3j. For the purpose of this analysis, mod- 
ifications of the 1 5o AW interaction were accomplished by 
adjusting the cr-meson coupling constant g^/^ Ji2| . In 
this way 1 Sq np (nn) interactions with different np (nn) 
scattering lengths were generated. 
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FIG. 2. Cross sections for 8 np = 28.0° (9 n i = 28.0°, <9 n2 = 
83.5°, 0i2 = 180°). The points are the data from the present 
work. The curves are MC simulations based on nd calcula- 
tions made with four values of a np : -22.0, -23.0, -23.75 and 
-25.0 fm. 



Simulated cross sections in comparison to our data for 
28.0° are shown in Figs. ^ and || for several values of a np 
and a nn , respectively. A value of a^N and its statistical 
uncertainty were determined for each detector-pair con- 
figuration using a single-parameter (either a np or a nn ) 
minimum-^ 2 fit to the absolute cross-section data. The 
results are given in Table [j]. The uncertainties listed in 
Table [| are statistical only. The systematic uncertainties 
in our determinations are ±0.8 fm for a np and ±0.6 fm 
for a nn . Uncertainties in the neutron detector efficien- 
cies and the integrated target-beam luminosity account 
for about 80% of the systematic uncertainty. 
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FIG. 3. Cross sections for 9 nn = 28.0° (0„i = 8 n2 = 28.0°, 
0i2 =0°). The points are the data from the present work. 
The curves are MC simulations based on nd calculations made 
with three values of a nn : -17.7, -18.0, and -19.0 fm. 

We observed no signficant angle dependence in a np , 
and the consistency in the a nn data from one angle to 
the next is statistically acceptable. Combining the sta- 
tistical and systematic uncertainties in quadrature, we 
obtain a np — -23.5 ± 0.8 fm and a nn — -18.7 ± 0.6 fm. 
Our result for a np is in agreement with the value of a np 
(-23.748±0.009 fm j|) obtained from free np scatter- 
ing measurements. We use this result to set an up- 
per limit on the influence of 3NF on the value of AW 
scattering lengths determined from our experiment, i.e., 
Aal% F < - a f np ee = 0.2 ± 0.8 fm, where free and nd 
refer to values obtained from data for free np scattering 
and from nd breakup, respectively. This result is consis- 
tent with zero. Scaling our result for a np by the ratio of 
a nn to a np , we obtain the upper limit due to 3NF effects 
in the nd breakup reaction to be Aa^ F < 0.2 ± 0.6 fm, 
which is also consistent with zero. 



Using the TM 3NF model we estimated possible effects 
of 3NF on np FSI cross sections. The calculations were 
produced by solving Eq.(|^) using four modern iVTV poten- 
tials: AV18 |§, CD Bonn @, Nijml and Nijmll Q. In 
each calculation the 27r-exchange TM 3NF potential |t]] 
was included as V4, which was split into 3 parts where 
each one was symmetrical under exchange of two par- 
ticles. For instance, for the 27r-exchange 3NF [[tJ, this 
corresponds to the three possible choices of the nucleon 
which undergoes the (off-shell) tt — N scattering. In our 
calculations with the 3NF the strong cut-off parameter A 
in the TM 3NF model was adjusted separately for each 
NN potential to reproduce the experimental triton bind- 
ing energy Jl(|. For details of the formalism and the 



numerical treatment refer to refs. 1Q 17]]. We calculated 
the percentage deviation caused by effects of the TM 3NF 
on the np FSI point-geometry cross section as the func- 
tion of 6 np at an incident neutron energy of 13 MeV. At 
all angles and for all potentials the change in the cross 
section due to the addition of the TM 3NF never exceeds 
6%. For the np FSI production angles of the present ex- 
periment it is limited to the range of 1% to 4%, which 
corresponds to a theoretical range of (Aa^ F )th from - 
0.8 to -0.2 fm. Our experimentally determined Aa^ F is 
within two standard deviations of the predictions using 
any of the four NN potentials with the TM 3NF adjusted 
to fit the triton binding energy. 



TABLE I. The a np and a nn values extracted from the fit to 
the present data for the absolute AW FSI cross section in nd 
breakup at the angles measured in the present study. The \ 2 
per datum for the best is given at each angle. The weighted 
mean of the data for all angles is given in the bottom row. 
All uncertainties are statistical only. 
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-18.7 ± 0.1 
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Summarizing, our measured values -23.5 
±0.8 fm and a nn — -18.7 ±0.6 fm. In our analysis mag- 
netic interactions are not considered. Their possible ef- 
fects have to be studied. By comparing our results 
for a np to the recommended value from np free scatter- 

nn to a np , we set an 



ing and scaling by the ratio of a, 

upper limit of Aa^ l F — 0.2 ± 0.6 fm on the contribution 
of 3NF effects on our value of a nn . Though the experi- 
mental results suggest an opposite sign for Aa'f l 1 p fF than 
predicted using the TM 3NF, our value is consistent with 
the theoretical predictions within the reported uncertain- 
ties. Since our value for a np obtained from nd breakup 
agrees with that from free np scattering, we conclude that 



our investigation of the nn FSI done under identical con- 
ditions should lead to a valid measure of a nn . Our value 
for a nn is in agreement with the recommended value [||, 
which comes from 7r~ d capture measurements, and dis- 
agrees with values obtained from earlier nd breakup stud- 
ies 01. 
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